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Nucleon transfer reactions in low-energy deuteron-deuteron scattering are described by solving 
exact four-particle equations in momentum space. The Coulomb interaction between the protons 
is included using the screening and renormalization method. Various realistic potentials are used 
between nucleon pairs. The energy dependence of the differential cross section, analyzing powers, 
polarizations, spin-transfer coefficient, and the quintet suppression factor is studied. 
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Neutron (n) and proton (p) transfer reactions in 
deuteron-deuteron (d-d) scattering, d + d —>■ p + ^ii and 
d + d ^ n + "^He, are amongst the simplest nuclear re- 
actions where charge-symmetry breaking (CSB) in the 
nuclear force can be verified or searched for. For that 
purpose one needs a very precise calculation of the four- 
nucleon (4A^) problem with different interaction models 
based on nucleon-nucleon {NN) and many-nucleon forces 
together with the inclusion of the Coulomb force between 
protons which is the most important cause of CSB. Such 
task is now possible in view of the progress achieved in the 
past few years [l|, 01 on the solution of^exact Alt, Grass- 
berger, and Sandhas (AGS) equations [3] for four-particle 
transitions operators that, in addition to the strong nu- 
clear force, include also the Coulomb interaction. 

The aim of the present paper is to study the energy 
dependence of the d + d —>■ p + ^ii and d + d ^>- n + '^He 
observables below three-body breakup threshold using 
different nuclear force models. In addition to the dif- 
ferential cross sections j^-Q there are precise measure- 
ments of the deuteron analyzing powers between 1.5 and 
4 MeV deuteron lab energy [y, Q . The analyzing powers 
show a complex structure of maxima and minima that in 
some cases varies rapidly with the energy. These rapid 
variations and complex structure means that different dd 
partial waves play a significant role, and agreement or 
disagreement with experimental data at different ener- 
gies may result from a delicate interplay between them. 
The contribution of each partial wave to the observables 
and their evolution with the energy is analyzed. 

Other observables that were measured but never ana- 
lyzed through exact 47V calculations are polarization of 
the outgoing nucleon 1 91 an d deuteron to nucleon spin 
transfer coefficients [10, Illl ; they are also calculated in 
the present work. 

Furthermore, we will get back to an old issue that re- 
mains of interest to d-d fusion in hot plasma. It was ini- 
tially believed T2] that at very low energy the production 
of neutrons resulting from d + d —>■ n + "^He could be con- 
trolled by polarizing all deuterons in the plasma such that 



only spin-aligned d-d fusion would take place. If spin- 



aligned d-\- d ^■ 



He reaction would be significantly 
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suppressed compared to d-\-d — >■ p-\-^il, one would have a 
significant reduction of undesired neutrons coming out of 
the plasma. In the past many approximate calculations 
were made that yielded contradicting results. More so- 
phisticated calculations jl3| indicated that no strong sup- 
pression should be expected. The present well-converged 
AN calculations aim at settling the assumption on the 
possible suppression of d + d — s> n + ■^He reaction with 
spin-aligned deuterons. 

Our description of the 4A^ system is based on exact 
four-particle equations for the transition operators as de- 
rived by Alt, Grassberger, and Sandhas [31; they are 
equivalent to the Faddeev-Yakubovsky equations lj| for 
the wave-function components. Since in the isospin for- 
malism protons and neutrons can be considered as identi- 
cal particles, the symmetrized form of the AGS equations 
[l|, 0] is appropriate. Although the initial d-d state has 
total isospin T = 0, the final n-'^He and p-^H states are 
dominated by both isospin 7" = and T = 1 compo- 
nents; a very small admixture of 7" = 2 is present due to 
the charge dependence of the hadronic and electromag- 
netic interactions. The most important source for the to- 
tal isospin nonconservation is the Coulomb force between 
the protons; it is included using the method of screening 
and renormalization |15l - [l8| . The hadronic charge depen- 
dence, as given by the modern NN potentials, is taken 
into account as well. After partial wave decomposition, 
the AGS equations are a system of three- variable integral 
equations that are solved numerically without any ap- 
proximation beyond the usual discretization of momen- 
tum variables on a finite mesh; technical details are given 
in Refs. [ll, Q. The results we present are well converged 
with respect to the partial-wave expansion, the number 
of momentum meshpoints, and the Coulomb screening 
radius used to calculate the short-range part of the am- 
plitudes. In this way the discrepancies with the experi- 
mental data may be attributed solely to the underlying 
NN forces or lack of many-nucleon forces. 

As two-nuclcon interactions we use the phenomeno- 
logical potentials Argonne V18 (AV18, hadronic part 
only) '19'], charge-dependent Bonn (CD Bonn) [2(|, and 



inside nonlocal outside Yukawa (INOY04) potential by 
Doleschall [2l|, and the one derived from the chiral per- 
turbation theory at next-to-next-to-next-to-leading order 
(N3L0) 22] . Furthermore, we consider also a two-baryon 
coupled-channel potential including virtual excitation of 
a nucleon to a A isobar [23] that in the AN system yields 
effective 3 A^ and AN forces [24!] . Point Coulomb is added 
for pp pairs. '^He and '^H binding energies calculated with 
those potentials are collected in the Table HI Thus, the 
presence of the 3N force is also simulated by the potential 
INOY04 that fits both '^He and '^H experimental binding 
energies and thereby of all the used potentials is the only 
one that reproduces correctly the momenta of the final 
n-^He and p-^H states. For this reason it is not surprising 
to see that INOY04 potential gives the best description 
of d+d -^ n + ^He and d+d ^ p + ^H data. We therefore 
show the predictions of all potentials for the reactions at 
deuteron lab energy Ed = i MeV but study the energy 
dependence of the observables using INOY04 only. 
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FIG. 1. (Color online) Differential cross section of d -I- d — >■ 
p -I- H and d + d ^ n + He reactions at 1.5, 3, and 4 MeV 
deuteron lab energy as function of the nucleon cm. scatter- 
ing angle. Results obtained with various realistic JVA'^ poten- 
tials are compared with the experimental data from Refs. [J] 
(squares), [a, 0] (circles), and [3 (triangles); the latter set is 
taken at Ea = 3.7 MeV. 



TABLE I. ^H and ^He binding energies (in MeV) and 
deuteron D-state probability Poid) (in percent) for different 
A^TV potentials. 

Although d-d elastic scattering is calculated simulta- 
neously with the transfer reactions, it is a less interesting 
case and therefore not discussed in the present work: As 
demonstrated in Refs. @,[2J]' the d-d elastic cross section 
data is well described by all NN force models and the 
corresponding deuteron analyzing powers are very small 
with quite large error bars that preclude physics conclu- 
sions. 

In Fig.[T]we present the differential cross section da/dil 
results for both d + d ~> p -{- ^H and d + d ^ n + ^He re- 
actions at Ed = 1.5, 3, and 4 MeV as function of the nu- 
cleon center-of-mass (cm.) scattering angle 9c. m.- This 
observable is symmetric with respect to 8c. m. — 90°. 
While at the lowest energy the differential cross section 
is flat between 0c. m. = 60° and 120°, with increasing 
energy it develops a minima around 0c. m. = 55° and 
125° and a local maximum at 0c. m. = 90°. Furthermore, 
da/dCl increases with energy at forward and backward 
angles. There is a good agreement between the experi- 
mental data and the predictions of the INOY04 potential 
whereas other models overestimate the data, especially at 
forward and backward angles; furthermore, the discrep- 
ancy seems to be proportional to the defect in theoreti- 
cal '^He and '^H binding energies. Thus, the differential 
cross section in the considered transfer reactions corre- 
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FIG. 2. (Color online) Deuteron vector analyzing power iTn 
of d + d — )■ p + '^H and d + d — > n-|-^He reactions at 1.5, 3, and 
4 MeV deuteron lab energy. The data are from Ref. [y, 01 for 
d + d -J- p -I- ^H and from Ref. [I ioT d + d^n + ^He. 



lates nearly linearly with 3N binding energies. 
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FIG. 3. (Color online) Same as Fig. [2] but for the deuteron 
tensor analyzing power T2o- 



FIG. 5. (Color online) Same as Fig. [2] but for the deuteron 
tensor analyzing power T22- 
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FIG. 4. (Color online) Same as Fig. [2] but for the deuteron 
tensor analyzing power T2i. 



In Figs. [2] - [5] we present the corresponding results 
for the deuteron analyzing powers. The vector analyz- 
ing power zTii in Fig. [2] varies quite slowly with the en- 
ergy. The experimental data are well reproduced by the 
INOY04 model, although a small discrepancy remains 
around 8c. m. = 45°. More significant discrepancies take 
place in the tensor analyzing power T2Q (Fig. [3]) at small 



scattering angles, although this observable is slightly 
overestimated almost in the whole kinematical regime. 
Nevertheless, the energy dependence of T20, being con- 
siderably stronger than for zTn, is well reproduced by the 
calculations. The tensor analyzing power T21 (Fig. U), 
though showing strong energy dependence at intermedi- 
ate angles, is quite well accounted for by the INOY04 pre- 
dictions. The description of the tensor analyzing power 
T22 (Fig.[S|) is even better; this observable, much like iTu, 
varies slowly with the energy. The results obtained with 
other NN potentials at ii^^ = 3 MeV deviate from the 
respective data more significantly. Furthermore, even if 
the reproduction of experimental binding energies is an 
important factor, it is certainly not the only one that 
matters since the linear correlation with 3iV binding en- 
ergies seems to be violated, e.g., for iTu at 6c. m. be- 
tween 40° and 90°, or for T22 between Gem. = 60° and 
120°, the predictions of N3L0, CD Bonn, and CD Bonn 
+ A almost coincide despite quite significant differences 
(up to 0.4 MeV) in the respective 3iV binding energies 
while the predictions of INOY04 and AV18 stay well sep- 
arated. A closer look into the calculated properties of 
deuteron, "^He, and ■^H as given in Table U suggests that 
the considered spin observables correlate, in addition to 
the 37V binding energy, also with the Z?-state probability 
of deuteron Pd (d) . The correlation between Pd {d) and 
3N binding energy takes place for most phenomenologi- 
cal A^A^ potentials, but N3L0 and CD Bonn + A models 
clearly violate that correlation thereby allowing to study 
the dependence of 4A^ observables on both 3A^ binding 
energy and Poid). The results in Figs. [2] - [5] indicate 
that increasing the ZN binding energy and decreasing 
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FIG. 6. (Color online) Outgoing nucleon polarization Py of 
d + d -7- p + '^H and d + d -^ n+ He reactions at 2,3, and 4 
MeV deuteron lab energy. The data are from Ref. [^. 



Poid) niove the theoretical predictions into the same di- 
rection but the corresponding rates, i.e., strength of the 
correlations, depend on the observable and kinematical 
regime. One may conjecture that by a slight decrease of 
Poid) of INOY04 while keeping the ^He and ^H binding 
energies unchanged one might be able to cure the zTn 
discrepancy around 6c. m. — 45°. However, T20 shows 
quite weak correlation with Poid) and therefore its de- 
scription would not be improved significantly by a small 
decrease of Pd (d) . 

In Fig. ini we present the results for the polarization 
Py of the outgoing nucleon in d + d — > p-l-'^H and 
d + d ^ n + ^Ke reactions; this observable is equivalent 
to the nucleon analyzing power Ay in the time-reverse 
reactions p + ^H — > d + d and n + '^He -^ d + d. Py is 
antisymmetric with respect to 0c. m. = 90°; we therefore 
show only the angular regime up to 0c. m. — 90° that 



contains all the available data. Py shows quite strong 
energy dependence; in the d + d ^^ p + ^H case it even 
changes the sign when Ed varies from 2 MeV to 4 MeV 
being almost zero at Ed — 3 MeV. In the d+d —> n-t-^He 
case the observable changes with the same trend but with 
about 1.5 MeV shift in the energy. The qualitative repro- 
duction of the experimental data having large errorbars 
is quite successful by the potential INOY04 whereas pre- 
dictions of the other interaction models are further away. 
As the deuteron analyzing powers, the nucleon polariza- 
tion correlates not only with the 3A^ binding energy but 
also with Poid) as one can see by comparing N3L0, CD 
Bonn, and CD Bonn + A results. 

Next we consider double polarization observables for 
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FIG. 7. (Color online) Deuteron to nucleon spin transfer co- 
efficient Ky for d + d ^>- p + '^'H and d + d -^ n + ''He reactions 
at Gem. ~ 0° as function of the deuteron lab energy. The 
data are from Ref. [ip] (circles) and from Ref. \V^ (squares); 
both sets refer to the d + d — >■ n + ^He reaction. 



which, unfortunately, the experimental data are much 
scarcer. The deuteron to neutron spin transfer coefficient 
Ry has been measured in several experiments [lOl 111] , 
however, only for the neutrons emitted in the forward 
direction 0c. m. = 0°; the data for the corresponding ob- 
servable in the d + d -^ p + ^H reaction is not avail- 
able below the three-body breakup threshold. In Fig. [7] 
we compare the deuteron to neutron spin transfer co- 
efficient Ky{0°) calculated using the INOY04 potential 
with two sets of experimental data in the energy range 
0.75 < Ed < 4 MeV. There is a good agreement above 
Ed = 2 MeV while at lower energies, where the observ- 
able shows stronger energy dependence, the data points 
from Ref. [ll| are slightly overpredicted. We include in 
Fig. [7] also the calculated Ky{0°) ior the d + d^p + ^H 
reaction that shows a similar behavior. The dependence 
on the NN force model is quite weak for i^^(0°), con- 
siderably smaller than the experimental error bars: at 
Ed = 3 MeV the predictions of AV18, N3L0, CD Bonn, 
CD Bonn + A, and INOY04 are 0.622, 0.626, 0.630, 
0.632, and 0.633, respectively. 

In Fig. \8\ we study the contribution of various initial 
and final state partial waves characterized by the relative 
two-cluster, i.e., d-d, n-^He, and p-^H, orbital angular 
momentum L. Since d + d ^ p + ^H and d + d ^ n + ^He 
observables show similar behavior, we restrict our study 
to the latter reaction at Ed = 1.5 and 4 MeV. We per- 
formed a series of calculations with INOY04 potential 
including L < L^i^y^ with Lmax ranging from 1 to 5. The 
predictions with Lmax = 4 and 5 are very close, indicat- 
ing that the L = 5 contribution is very small; however, 
since the convergence is not monotonic, we have also ver- 
ified that L — 6 contribution can be safely neglected. In 
contrast, the partial waves with L — 4 yield a sizable 
contribution a.t Ed = 4: MeV and for the deuteron ten- 
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FIG. 8. (Color online) Observables of d + d ^ n + He re- 
action at Ed = 1.5 and 4 MeV. Results including initial and 
final states with two-cluster relative orbital angular momen- 
tum L < Lmax are compared for Lmax ranging from 1 to 5. 
imax = 4 and 5 results lie almost on top of each other. The 
experimental data are as in previous figures. 



sor analyzing powers even at Ed 
observables at this energy Lmax 



1.5 MeV. For other 
3 is sufficient while 



Py is reasonably well converged with L„ 



2. It is in- 



teresting to note that L,„ax = 1 results for most observ- 
ables vary slowly with the energy while stronger energy 
dependence comes from L — 2, 3, and 4 partial waves. 
The above analysis partially explains sizable differences 
between our converged results and those of Ref. '23 ob- 
tained using the resonanting group method (RGM) since 
the latter work included only L < 3 states. Our AV18 re- 
sults calculated with L,nax = 3 (not shown here) are qual- 
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FIG. 9. (Color online) The quintet suppression factor 0-2/0- 
for d + d -i- p + ^H and d + d — >■ n + '^He reactions as function 
of the deuteron lab energy. 



itatively similar to the corresponding results of Ref. [231 ■ 
Finally we present results for the observables char- 
acterizing the spin correlations of the initial-state 
deuterons. We are not aware of the existence of the ex- 
perimental data, but there are plans [2a, [23| for future 
experiments involving the direct measurement of the so- 
called quintet suppression factor (QSF) defined as the 
ratio 0'2/cr, where a2 is the cross section for the consid- 
ered transfer reaction with the spins of the initial state 
deuterons being parallel, i.e., with the total spin being 
2, and a is the unpolarized (spin- averaged) cross section. 
As mentioned in the introduction, this observable is rele- 
vant for d-d fusion in hot plasma. Results obtained using 
INOY04 potential for Ed between 50 keV and 4 MeV are 
shown in Fig. O The QSF for the d + d ^>- n + ^Re reac- 
tion is indeed smaller than for d + d — >p + ''H; however, 
the difference is only about 20 - 25% in the whole con- 
sidered regime. The energy dependence of the QSF is 
in both reactions similar: it is of the order of 0.25 above 
Ed = 2 MeV but increases more rapidly as the energy de- 
creases. Nevertheless, our QSF prediction are somehow 
smaller than the values obtained from i?-matrix analysis 
and RGM [l^l . This observable shows also strong depen- 
dence on the NN force model, e.g., the QSF values for 



the d-\- d 



'He reaction dX Ed = i MeV predicted 



by AV18, N3L0, CD Bonn, CD Bonn + A, and INOY04 
models are 0.133, 0.159, 0.164, 0.185, and 0.219, respec- 
tively. Thus, the QSF correlates with both 3iV binding 
energy and Pd (d) , increasing when the former increases 
and/or the latter decreases. 

In summary, we performed exact four-particle calcu- 
lations oi d + d ^i- p + ^H and d + d ^ n + ^He reac- 
tions with several realistic NN potentials. Energy de- 
pendence of the differential cross section and spin ob- 
servables was studied below three-body breakup thresh- 
old. Correlations of the predictions with 3iV binding en- 
ergy and deuteron Z?-state probability Poid) were ob- 



served. The INOY04 potential that fits both ^He and 
■^H experimental binding energies and has the smallest 
Poid) = 3.60% among all realistic potentials, accounts 
well for the experimental data with only few discrepan- 
cies, e.g., the one in the deuteron tensor analyzing power 



T20. Some data suggest that even smaller Poid) would 
be preferred. We also predict the quintet suppression 
factor in the d + d ^ n + '^He reaction to be only up to 
25% smaller than the one in d + d — s- p + "^H; the d-d fu- 
sion with spin- aligned deuterons seems to be significantly 
suppressed at few MeV but not in the keV regime. 
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